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16  Abstract 

Improved  azimuthal  resolution  of  proximate  aircraft  necessary  to  support  ATC 
automation  can  be  achieved  by  beacon  surveillance  systems  employing  monopulse 
angle  estimation  techniques  described  in  this  report . 

Included  in  the  report  are  the  results  of  beacon  surveillance  monopulse 
system  analyses  relating  to  of f-boresight  angle  estimation  using  short  (1/2  psec) 
pulses;  the  effects  of  specular  and  diffuse  multipath  signal  return;  the  effects 
of  overlapping  ATCRBS  fruit  replies,  and  the  problems  of  antenna  pattern  design. 
These  topics  have  been  studied  in  detail  as  part  of  the  Lincoln  Laboratory 
design  of  the  Discrete  Address  Beacon  System  (DABS) . 

This  report  summarizes  significant  analytical  results  obtained.  In  general, 
it  has  been  concluded  that  the  ATC  environment  does  not  pose  a serious  problem 
to  the  use  of  the  monopulse  concept  for  beacon  system  direction  finding  and  that 
sufficient  direction  finding  accuracy  can  be  obtained  using  a small  number  of 
narrow  pulses  for  each  scan. 
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DABS  MONOPULSE  SUMMARY  REPORT 


1.0  INTRODUCTION 

Studies  have  shown  [1]  that  beacon  surveillance  performance  required  to 
support  ATC  automation  can  be  achieved  by  minimizing  self  interference  and  by 
improving  resolution  of  proximate  aircraft.  Interference  can  be  reduced  by 
lowering  the  interrogation  rate.  Froximate  targets  can  be  resolved  by  interro- 
gating each  target  separately  (by  an  addressed  interrogation)  so  as  to  separate 
their  replies  in  time.  There  is  sufficient  time  on  the  up  and  down  links  to 
discretely  address  all  targets,  even  in  high  traffic  densities,  if  position  and 
identity  can  be  derived  from  about  two  replies  per  scan. 

Azimuth  estimation  in  the  current  ATCRBS  is  based  upon  a beam-splitting 
or  sliding  window  process.  The  estimate  accuracy  for  these  techniques  is 
related  to  the  beanwidth/runlength  ratio  using  runlengths  of  16  to  20  per  scan. 
A method  that  measures  azimuth  accurately  on  the  basis  of  a few  replies  per 
scan  is  required. 

This  capability  can  be  achieved  using  a monopulse  angle  estimation  system. 
Important  monopulse  considerations  are:  of f-boresight  angle  estimation  using 
short  (1/2  ysec)  pulses;  the  effects  of  specular  and  diffuse  multipath  signal 
return;  the  effects  of  overlapping  ATCRBS  fruit  replies,  and  the  problems  of 
antenna  pattern  design.  These  topics  have  been  studied  in  detail  as  part  of 
the  Lincoln  Laboratory  design  of  the  Discrete  Address  Beacon  System  (DABS). 

This  report  summarizes  the  significant  analytical  results  obtained.  In  gene- 
ral, it  has  been  concluded  that  the  ATC  environment  does  not  pose  a serious 
problem  to  the  use  of  the  monopulse  concept  for  beacon  system  direction  finding 
and  that  sufficient  direction  finding  accuracy  can  be  obtained  using  a small 
number  of  narrow  pulses  for  each  scan. 

Sections  2.,  3.,  and  4.  of  the  report  describe  antenna  characteristics 
required  to  support  monopulse,  and  summarize  the  features  of  several-- possible 
hardware  realizations  of  the  monopulse  processor.  The  realization  selected  is 
the  so-called  half-angle  phase  comparator  which  utilizes  phase  detectors  to 
produce  an  unambiguous  of f-boresight  indication  over  a large  fraction  of  the 
beamwidth.  Sections  3.  and  4.  discuss  the  effects  on  monopulse  performance  of 
interference  from  various  sources  and  of  various  antenna  design  parameters. 

The  impact  of  diffraction  of  the  incoming  wavefront  around  obstructions  such 
as  nearby  buildings  is  also  briefly  discussed. 

Sections  5.  through  7.  describe  the  performance  actually  obtained  using 
the  monopulse  processor  employed  at  DABSEF , an  implementation  of  the  processor 
specified  in  the  DABS  sensor  engineering  requirement.  Section  5.  focuses  on 
the  stability  of  the  processor  and  its  sensitivity  to  the  frequency  and  ampli- 
tude of  the  received  signal.  Section  6.  describes  the  method  used  to  evaluate 
reply  accuracy  as  a function  of  of f— boresight  angle,  and  gives  results  for  a 
large  number  of  ATCRBS  targets  of  opportunity.  In  Section  7.  target  report 
accuracy  is  evaluated  using  smoothing  techniques  to  estimate  aircraft  position. 
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Section  8.  describes  analytical  and  experimental  work  done  at  Lincoln  Laborato 
to  quantify  and  confirm  the  effects  of  obstacle  shadowing  upon  DABSEF  mono- 
pulse direction  finding  accuracy.  Conclusions  are  drawn  in  Section  9. 

Appendix  A describes  the  system  used  to  evaluate  DABS  monopulse  perform- 
ance (DABSEF) . An  error  analysis  is  included  to  estimate  the  effect  on 
performance  of  variations  in  half-angle  processor  hardware  elements. 

Appendix  B is  a summary  of  a curve  fitting  technique  using  look-back 
and  look-ahead  data  to  accurately  develop  aircraft  trajectory.  This  tech- 
nique is  the  basis  for  calculating  the  DABSEF  monopulse  angle  estimator's 
accuracy . 
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2.0  MONOPULSE  PROCESSING 


2.1  Introduction 

The  ability  of  monopulse  techniques  to  estimate  target  azimuth  from  a 
single  received  pulse  is  based  on: 

a.  An  antenna  that  can  provide  outputs  related  to  the  angle  of  in- 
cidence of  the  incoming  plane  wave  (of fboresight  angle), 

b.  A processor  that  can  convert  the  antenna  outputs  to  a signal  re- 
lated to  of fboresight  angle,  and 

c.  A calibration  method  which  can  establish  the  exact  relation  be- 
tween processor  output  and  actual  offboresight  angle. 

Section  2.2  introduces  the  subject  of  monpulse  antennas,  and  notes  their 
principal  parameters  and  requirements.  Detailed  discussion  of  monopulse  an- 
tennas, their  influence  on  overall  monopulse  angle  estimation  accuracy,  and 
characteristics  which  minimize  their  vulnerability  to  interference  is,  how- 
ever, deferred  to  Section  4.  Section  2.3  proceeds  with  a detailed  discu- 
sion  of  the  monopulse  processor,  noting  differences  between  amplitude  and 
phase  comparison  processing  and  explaining  why  the  so-called  half-angle  phase 
comparator  has  been  selected  for  the  DABS  application.  For  purposes  of  this 
discussion,  the  monopulse  antenna  is  treated  as  a simple  "three-port  network" 
providing  sum  (£),  difference  (A)  and  omnidirectional  (ft)  control  outputs. 

2.2  Monopulse  Antennas 

A monopulse  antenna  requires: 

- A "Sum"  pattern  (£)  corresponding  to  a symmetric  directional  mainlobe, 
typically  a few  degrees  wide. 

- An  asymmetric  "Monopulse  Difference"  pattern  (A)  with  a directional  pat- 
tern commensurate  in  width  with  the  £ pattern  and  accurately  centered  with  re- 
spect to  it.  Its  signals  are  used,  in  conjunction  with  those  from  I,  to  determine 
the  bearing  angle  of  targets  known  to  be  in  the  main  beam  (sometimes  referred 

to  as  "monopulse  window"). 

- A "Control"  pattern  ( ft,  often  implemented  as  an  omnidirection  pattern). 
This  pattern  is  used  in  conjunction  with  the  £ pattern  to  provide  the  various 
transmit  sidelobe  suppression  functions  (SLS)  and  the  receive  sidelobe  flagging 
functions  (RSLS) . 

Fig.  2-1  is  a simple  representation  of  a monopulse  antenna  in  which  the 
antenna  outputs  correspond  to  these  three  essential  patterns. 

The  types  of  antennan  under  discussion  when  operating  in  conjunction  with 
the  monopulse  processors  discussed  in  the  following  section  introduce  system 
angle  errors  which  are  dependent  on  the  beamwidth  and  on  the  difference  pattern 
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Fig.  2-1.  Monopulse  antenna  (simplified  representation) 


shape.  It  is  common  practice  to  characterize  the  slope  of  the  difference  pattern 
by  the  location  of  the  sum-to-difference  cross-over  point  (A/Z  = 1);  this 
essentially  specifies  the  difference  pattern  beamwidth  relative  to  that  of  the 
sum.  The  desirable  location  of  the  cross-over  point  varies  between  -3  dB  on 
the  sum  pattern  (this  makes  the  aperture  width  required  for  the  difference 
compatible  with  that  required  for  the  sum)  and  -4  dB  (anything  less  implies 
that  the  aperture  is  inefficiently  utilized) . This  desirable  range  of  cross- 
over values  results  in  a small  variation  in  accuracy;  the  sum  beamwidth  emerges, 
therefore,  as  the  dominant  antenna  determinant  of  inherent  accuracy.  Lincoln 
Laboratory's  experience  with  the  monopulse  receiver  and  antenna  at  DABSEF 
indicates  that  it  is  reasonable  to  expect  field  deployable  equipment  to  achieve 
a beamsplit  factor  of  40:1.  Thus,  a 4 sum  azimuth  beamwidth  is  about  the 
upper  limit  if  inherent  direction  finding  accuracy  of  0.1  rms  is  desired. 

The  relative  importance  of  monopulse  antenna  design  parameters  upon  mono- 
pulse processing  performance  and  upon  overall  monopulse  angle  estimation  system 
performance  is  discussed  in  detail  in  Section  4. 

2.3  Monopulse  Processors 

The  monopulse  processor  provides  a signal  indicative  of  the  azimuth  angle 
between  the  target  and  the  antenna  pointing  direction  for  each  received  pulse. 

Two  basic  receiver  configurations  can  be  used  to  generate  "monopulse"  signals: 
amplitude  comparison  and  phase  comparison  [3]. 


In  the  amplitude  comparison  scheme,  shown  in  Fig.  2-2a  (sometimes  refer- 
red to  as  the  Chubb  approach)  RF  antenna  outputs  corresponding  to  two  angle- 
squinted  beams  are  logarithmically  amplified  and  detected,  and  then  subtracted, 
yielding  a bi-polar  video  from  which  of f-boresight  angle  can  be  deduced.  The 
angular  region  over  which  this  scheme  generates  an  unambiguous  output  is  more 
limited  than  desired.  This  is  true  since  when  the  squinted  beams  are  generated: 

(i)  directly  by  the  antenna,  the  angular  limit  is  determined  by  the 
first  nulls  of  the  individual  patterns. 

(ii)  from  linear  combinations  of  independent  sum  and  difference  beams 
(Z  + A,  and  Z - A) , the  limit  corresponds  to  the  sum  difference 
crossover  or  about  the  3 dB  beamwidth. 

This  limitation  reduces  the  desired  flexibility  for  DABS  interrogation  sche- 
duling. Specifically,  it  may  preclude  direction  finding  in  situations  when 
there  is  sufficient  signal  strength  to  perform  detection  and  communication, 
for  example,  for  near-in  targets  outside  the  3 dB  beamwidth. 

The  phase  comparison  scheme  shown  in  Fig.  2-2b,  (sometimes  referred  to  as 
the  Bell  Labs  approach) , had  its  origin  in  the  simple  direction  finding  scheme 
in  which  the  bearing  angle  is  obtained  by  measuring  the  relative  phase  between 
two  displaced  antennas.  This  scheme  is  refined  in  optimized  monopulse  systems  by 
first  generating  independent  sum  and  difference  patterns,  each  having  low  sidelobes. 
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and  making  maximum  use  of  the  available  aperture,  then  combining  the  sum  and 
difference  outputs  in  a hybrid  to  yield  outputs  Z + jA  and  £ - jA.  It  can  be 
verified  that  the  aperture  illuminations  associated  with  these  new  patterns 
tend  to  look  like  overlapping  antennas  with  displaced  phase  centers.  The 
normalization  of  the  two  signals  is  accomplished  by  phase-matched  limiters. 

The  bi-polar  video  output  of  the  phase  detector  contains  all  the  angle  inform- 
ation. 

It  has  been  shown  in  [4]  that  the  optimum  azimuth  estimate  in  a receiver 
noise  background  can  be  obtained  by  solving 


V6> 

gec§> 


= tan 


_ sin  <t> 

2 1 + cos  <(> 


(2-1) 


<f>  = arg(E  + jA)  - arg(E  - jA) 


(2-2a) 


0 = offboresight  angle  estimate 


(2-2b) 


and  where  Gj;(0) , G^(0)  represent  the  amplitudes  of  the  sum  (E)  and  difference 
(A)  radiation  patterns.  The  ratio  G^(0)/Gj;(0)  is  often  referred  to  as  the 
monopulse  function.  For  large  signal- to-noise  ratio  (SNR)  and  targets  near 
boresight,  it  turns  out  that|A/E|  <<  1 and  i}  is  a small  angle.  In  this  case 

cos  <p  ~ 1 so  that 


GA(6) 

ge(§) 


~ 1/2  sin  <J> 


(2-3) 


Therefore  the  detector  shown  in  Figure  2-2b  with  output  represented  by  a sim- 
ple sine  function,  represents  an  optimum  configuration  under  the  above  condi- 
tions. Its  output  can  be  fed  into  a table  look-up  to  invert  the  monopulse 
function  and  hence  provide  an  implementation  of  (2-3) . Unfortunately  as 
| A/E  J -+■  1 serious  noise  errors  begin  to  occur  and  then  for  | A/ E | > 1 the 
sinusoidal  phase  detector  output  becomes  ambiguous.  This  results  in  the  same 
restriction  in  the  unambiguous  monopulse  "field  of  view"  as  the  previous  am- 
plitude comparison  system. 


Fortunately,  unlike  the  amplitude  system,  the  phase  comparison  system 
can  be  made  to  work  over  the  full  width  of  the  beam  (where  ] A/ E | is  monotoni- 
cally  increasing).  Two  classes  of  design  options  are  available.  In  the  first, 
a quadrature  channel  is  provided  to  measure  cos  [arg(E  + jA)  - arg(E  - jA)]. 
This  can  be  used  by 
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(1)  taking  the  sine  output  as  the  primary  estimate  and  use  using  the 
polarity  of  the  cosine  output  to  resolve  the  sine  ambiguity. 

(2)  taking  the  output  of  the  phase  detector  that  has  the  more  favor- 
able characteristics  for  the  specific  measurement,  e.g.,  use  the 
sine  output  for  |sin  <j)|  <_  0.707,  cos  <j>  < 0,  and  use  the  cosine 
output  for  | cos  <j>|  < 0.707. 

(3)  using  both  outputs  together  to  compute,  or  look  up  in  a table, 
the  estimate  of  of f-boresight  angle  using  the  optimum  estimation 
equation  (2-1) . 

None  of  the  above  alternatives  is  entirely  satisfactory.  The  first  is  sub- 
ject to  large  errors  in  the  estimate  near  |a/Z|  = 1,  while  the  other  two  lead 
to  rather  complex  implementations. 

The  second  class  of  design  options  to  overcome  the  ambiguity  problem 
uses  a "half- angle  processor".  If  angles 


a,  6 = arg  (Z  ± jA)  - arg  Z 


(2-4) 


are  defined  and  it  is  noted  from  (2-2)  that  4>  = (a-  6),  (2-1)  becomes 


Gz(6) 


tan 


(a  - B) 
2 


(2-5) 


Since  the  angles  are  always  less  than  90°,  phase  detectors  are  available 

that  will  linearly  and  unambiguously  measure  the  phase  difference  between  the 
signals  £ ± jA  and  I.  Therefore  summing  the  outputs  of  the  two  half-angle 
phase  detectors  yields  an  azimuth  estimate  that  is  identical  to  the  optimum 
Bell  Labs  processor  without  any  inherent  ambiguity. 

A realization  of  the  half-angle  processor  is  shown  in  Figure  2-3.  Al- 
though the  two  pairs  of  limiter  channels  do  not  carry  equal  level  signals, 
the  amplitude  difference  will  not  be  large  because  of  the  quadrature  rela- 
tionship between  Z and  jA  and  the  fact  that  since  detectability  must  always 
be  maintained,  the  Z signal  cannot  vanish.  It  is  desirable  to  use  two  redun- 
dant half-angle  processors,  rather  than  one,  because  the  sensitivity  to  phase 
errors  between  the  channels  Z + jA  and  A + jZ  is  half  as^great  for  the  con- 
figurations shown  as  it  would  be  for  the  corresponding  channels  Z + jA  and  Z. 
Furthermore,  the  additional  phase  detector  yields  a 3 dB  improvement  in  SNR. 

As  suggested  by  (2-5)  it  remains  to  take  the  processor  output  (a-  B) 
and  calculate  tan  [ (ct  - B)/2],  which  is  used  with  an  inversion  of  the  mono- 
pulse curve  to  obtain  the  azimuth  estimate.  This  procedure  requires  separate 
calibration  of  the  monopulse  function  and  the  phase  detector  characteristics. 
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In  practice  it  is  easier  to  obtain  a calibration  curve  that  relates  the  azimuth 
estimate  directly  to  the  phase  detector  output  (a  - 3)  since  this  also  ac- 
counts for  receiver  nonlinearities.  Hence,  the  final  form  for  the  azimuth  es- 
timator is 


0 = f(a  - 6) 


(2-6) 


where  the  function  table  f(.)  is  generated  from  the  quantized  version  of 
(a  - 6) . 

To  complete  the  description,  Figure  2-4  shows  the  output  of  the  half-angle 
processor  vs.  A/E  for  two  types  of  phase  detectors.  The  relation  2 tan--*-  A/E 
pertains  to  an  ideal  triangular  or  sawtooth  phase  detector,  and 


2Real  A/E 


applies  to  a sinusoidal  device.  Practical  phase  detector  characteristics  would 
lie  between  the  two  curves.  The  actual  behavior  of  the  processor  output  vs. 
target  position  of f-boresight  is  also  a function  of  the  antenna  monopulse 
characteristic,  i.e.  , A/E  vs.  azimuth. 

In  summary,  an  azimuth  estimate  is  obtained  by  processing  the  A and E sig- 
nals with  the  monopulse-processor  circuit."*-  The  output  of  this  device,  refer- 
red to  as  monopulse  video,  MV,  is  a monotonic  function  of  the  complex 

J A 
ratio  j 


MV  = f (A/E) 


(2-7) 


For  the  case  when  there  is  no  noise  and  no  interference, 

| = E(0)  (2-8) 

where  E(9)  is  the  normalized  difference  pattern  as  a function  of  offboresight 
angle  0 


E (0) 


V6> 

Ve> 


and  E are  the 


difference  and  sum  beam  outputs. 


(2-9) 


Typical  G^(0)  and  G^(0)  characteristics  are  shown  in  Figure  2-5.  The 
magnitude  of  the  normalized  difference  pattern  of  the  DABSEF  array  antenna 
system  is  shown  in  Figure  2-6.  Note  (from  Figure  2-5)  that  E(0)  is  « real 
for  | 0 | < 3.0°  and  is  monotonically  increasing  for  this  range  of  0.  Thus, 
there  is  an  unambiguous  relationship  between  MV  and  0;  i.e., 


MV  = f [E (0) ] = h(0) 


(2-10) 


This  function  h(0)  is  referred  to  as  the  monopulse  calibration  curve+.  Be- 
cause h(0)  increases  monotonically,  it  possesses  an  inverse,  defined  as 

0 = h_1  (MV)  = g(MV)  (2-11) 


Thus,  the  azimuth  estimator  is  given  by  the  relation 


0 = g(MV)  = g[f (A/Z) ] 


(2-12) 


Equations  2.7-2.12  are  shown  graphically  in  Figure  2-7.  The  operational  pro- 
cedure for  obtaining  g(MV) , the  monopulse  calibration  table,  is  described  in 
Section  3.3  of  Appendix  A. 

2.4  Monopulse  Direction  Finding  Accuracy 

The  accuracy  with  which  the  azimuth  angle  of  a target  can  be  estimated  us- 
ing th'>  processors  described  in  Section  2.3  in  the  absence  of  multipath  and 
interference,  has  been  called  "inherent"  accuracy.  The  approach  that  has  been 
used  is  to  evaluate,  for  various  sources,  the  error  in  reported  azimuth  when 
only  a single  reply  is  available.  This  will  be  taken  as  the  "baseline  perform- 
ance." Rather  than  following  the  convention  of  resolving  overall  errors  into 
contributions  by  subsystems,  the  errors  will  be  ordered  according  to  their  sta- 
tistical nature. 


2.4.1  Totally  Random  Errors 


In  practice,  receiver  noise  is  the  only  source  of  error  which  is  random 
on  a pulse-to-pulse  basis.  The  standard  deviation  of  the  error  (a)  is  given 
by  the  following  formula  [5] 


1 

a = 

>'2  • N • SNR 


vC 


jdcel 


E’(0) 


(2-13) 


+This  relationship  may  be  determined  by  interrogating  a fixed  transponder  at 
an  accurately  surveyed  position.  See  Appendix  A,  Section  3.3. 


where 


SNR  = Signal  to  noise  ratio  (at  electrical  boresight) 
a = Standard  deviation  of  azimuth  error,  degrees 
0g  = Beamwidth  (3  dB) , degrees 
N = Number  of  independent  samples 
0 = Of f-boresight  angle  in  degrees 
E'(0)  = Slope  of  normalized  difference  pattern 

An  approximate  formula  assumes  a linear  E(0)  and  cross-over  at  the  3 dB 
point,  in  which  case  E(0)  = 20/0g.  It  is  then  appropriate  to  define  the 
quantity  0g/a,  which  is  sometimes  referred  to  as  the  "beam-split  ratio"  and 
is  plotted  in  Figure  2-8  as  a function  of  offboresight  angle  for  several  values 
of  the  product  SNR  x N (effective  SNR). 

In  cases  where  E(0)  is  non-linear,  noise  effects  introduce  a bias  error 
that  increases  in  magnitude  with  of f-boresight  angle.  Figure  2-9  shows  an 
analytical  result  for  bias  and  sigma  for  an  E(0)  typical  of  the  DABSEF  array 
antenna.  Note  that  the  noise  induced  bias  is  only  a fraction  of  the  noise  in- 
duced standard  deviation  (sigma). 

2.4.2  "Reply-Dependent"  Errors 

Certain  errors  are  constant  from  pulse  to  pulse,  within  a reply,  but 
their  magnitude  depends  on  some  of  the  parameters  which  describe  individual 
interrogations.  For  example,  error  magnitude  is  dependent  on  the  off- 
boresight angle  which,  for  a given  sweep,  can  be  thought  of  as  partially  ran- 
dom. However,  depending  on  the  design  of  the  interrogation  scheduler,  the 
offboresight  angle  may  not  be  truly  random.  When  several  replies  are  received 
on  a single  scan  at  different  offboresight  angles  one  can,  through  an  averag- 
ing process,  smooth  out  the  azimuth  estimate  in  a manner  analogous  to  the  way 
that  independent  pulse  samples  are  averaged  within  a reply  to  reduce  the  ef- 
fects of  noise. 

Calibration  errors  are  important  contributors  to  this  type  of  inaccuracy. 
When  a lookup-table  is  used  to  relate  processor  output  to  of f-boresight  angle, 
then  the  error  becomes, 


where 


AOtrue  = true  boresight  angle 

C (.)  = functional  relation  between  A0  true  and  processor  output 
for  a given  target  and  target  position 

C(.)  = stored  look-up  table 


feflJHKr 


1000 


i 
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BIAS  (DEG) 
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Insight  into  this  behavior  can  be  gained  by  assuming  linear  C(.)  and  CT(.) 


CT(y)  = ay 


C(y)  = a'  y + b 


, true  , , 

e = a' + b = A0^ 

a true 


,(-r-i)+b 

/a'  \ 

ds  onl lldue  tc 

cent  'a  ana  0.1°, 


Typical  upper  bounds  onl lldue  to  variations  in  signal  charac- 

teristics are  a few  percent  ' ancr  0.1°,  respectively.  Experimental  measure- 
ments of  ( &L  -A  for  many  beacon  targets  of  opportunity  appear  in  Appendix  A. 

\ 3 ” / 


The  effect  of  a difference  in  target  elevation  angle  (a),  compared  to  cali- 
bration source  elevation  angle  (ac)  can  be  derived  for  the  linear  region  of 
the  processor  from  geometric  considerations  with  the  result  that: 


a'  y cos  a 
a cos  a 


Figure  2-10  shows 


vs 


a for  various  values  of  a . 

c 


Three  significant  observations  may  be  made  with  respect  to  elevation  in- 
duced errors.  First,  the  "reply  bias"  errors  tend  to  be  of  opposite  polarity 
on  either  side  of  boresight;  this  provides  a rapid  reduction  of  the  net  error 
when  two  or  more  replies,  sufficiently  separated,  are  available.  Secondly, 
the  errors  are  essentially  independent  of  beanwidth;  the  driving  parameter  is 
of f-boresight  angle.  Thirdly,  although  azimuth  accuracy  is  usually  used  to 
characterize  position  measurement,  the  associated  cross-range  accuracy  is  for 
most  purposes  more  directly  pertinent.  Figure  2-11  shows  the  cross-range  error 
as  a function  of  ground  range  for  various  aircraft  altitudes.  The  curves  shown 
extend  inward  only  to  a ground  range  corresponding  to  about  40°  in  elevation 
and  indicate  that  the  cross-range  errors  vary  from  tens  to  a few  hundreds  of 
feet. 


Other  contributions  to  the  reply-dependent  errors  are: 

A/D  quantization  of  the  monopulse  signal:  for  an  8-bit  (7  plus  sign)  conver- 

ter, the  rms  error  is  about  1/4%  of  the  maximum  of f-boresight  angle.  Azimuth 
shaft  encoder:  typical  azimuth  shaft  encoder  errors  of  less  than  .02  degrees 

rms  are  achieved. 
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Fig.  2-11.  Calibration  induced  cross-range  error. 
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2.4.3  Scan-Dependent  Errors 

There  are  errors  which  vary  only  from  scan  to  scan  and  therefore  cannot 
be  eliminated  by  any  form  of  processing.  Wind-induced  deflections  of  the 
antenna  aperture  causing  an  electrical  boresight  shift  is  one  example.  It  has 
been  estimated  from  past  experience  that  the  maximum  shift  for  a 40  ft.  aper- 
ture can  be  kept  to  less  than  0.05o  within  the  usual  range  of  FAA  environmental 
conditions.  Another  such  error  can  be  caused  by  drive  gear  backlash  in  the 
presence  of  wind.  For  typical  drive  mechanisms,  this  error  has  been  estimated 
at  0.02°  maximum. 

Whereas  the  above  errors  are  related  to  time-varying  factors,  there  are 
scan  dependent  errors  which  are  also  azimuth-dependent.  For  example,  a ra- 
dome  can,  in  principle,  be  a source  of  such  an  error;  however,  the  use  of  a 
foam  radome  has  been  found  to  virtually  eliminate  this  problem.  The  rotary 
joint  can  also  introduce  an  error  if  the  £ and  A channels  do  not  track  in 
amplitude.  A 0.05  dB  differential  error  appears  to  be  a reasonable  design 
specification  and  leads  to  an  error  which  is  at  most  1/2%  of  the  of f-boresight 
angle. 


2.4.4  Scan-Independent  Errors 

The  errors  considered  here  are  system  biases.  Sources  of  such  errors 
include  errors  in  the  initial  system  boresight  alignment  during  calibration. 
None  of  these  appear  sufficiently  large  to  be  taken  into  account  in  the  total 
error  budget.  Also,  since  there  are  variations  in  the  antenna  boresight  as  a 
function  of  elevation  which  cannot  be  taken  into  account  in  the  calibration  or 
alignment,  this  causes  an  additional  bias  error  for  a target  which  does  not 
change  altitude  rapidly  from  scan-to-scan.  This  error  is  primarily  attributed 
to  aperture  errors,  and  is  proportional  to  beanwidth.  The  nominal  figure  is 
±0.01  per  degree  of  beamwidth,  maximum. 

2.4.5  Direction  Finding  Summary 

An  overall  summary  of  monopulse  direction  finding  accuracy  is  given  in 
Table  2-1  for  beamwidths  of  4°  and  2°.  The  numbers  tabulated  are  to  be  inter- 
preted as  achievable  performance.  Contributions  which  are  negligible  have 
been  omitted. 
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DIRECTION  FINDING  ACCURACY  SUMMARY  FOR  ROTATORS 
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3.0  EFFECTS  OF  INTERFERENCE  ON  MONOPULSE  PERFORMANCE 


3.1  Introduction 


In  addition  to  the  hardware  and  receiver  noise  errors  discussed  in  the 
preceding  section,  the  performance  of  the  monopulse  processor  is  affected  by 
the  presence  of  other  time  coincident  signals  such  as  ATCRBS  fruit  and  ground 
reflection  multipath.  Detailed  studies  of  these  effects  have  shown  that,  in 
addition  to  causing  an  increase  in  the  estimate  variance,  the  presence  of 
interference  can  cause  an  estimator  bias.  For  mainbeam  interference  [6]  it 
has  been  shown  in  that  this  bias  is 


— — p (p  + cos  <}>) 

(0  - e ) = (e  - e ) — 0 - 2 

1 + 2p  cos  4>  + p 
o 


(3-1) 


where 


h W 

P°  “ *s  W 


1 

SIR 


(3-2) 


where  A^,  0^  are  the  amplitude  and  azimuth  of  the  signal,  A^. , 0^.,  the  amplitude 
and  azimuth  of  the  interferer  and  <f>  is  the  relative  phase  between  the  DABS  and 
the  interference,  and  SIR  is  the  signal  to  interference  ratio.  Eqn  • 3.1  assumes 
the  half-angle  processor  (Eqn.  2.13)  is  linear  and  modelled  by 

f(£)  = Re(|)  (3-3) 

and  the  normalized  ratio  pattern  and  calibration  curves  are  also  linear. 

The’  bias  error  is  plotted  in  Fig.  3-1  for  a target  on  boresight,  and  an 
interferer  at  the'’T,  dB  beam  edge.  In  general,  the  curve  demonstrates  the 
large  bias  that  can  occur,  especially  at  the  out-of-phase  condition. 

One  can  more  accurately  model  the  half-angle  processor  as  non-linear 
by  Eqn.  2.13, 


f(f)  = Re  (|)| 


1 + lil2 


2 Im  (Y) 


1 + |£|2  + 2 Im  (|) 


(3. A) 


and  the  normalized  ratio  pattern,  E(0),  and  calibration  curve  as  non-linear. 


Figure  3-2  shows  the  difference  between  the  linear  and  the  non-linear 
approximations  for  pq  = -3  dB,  0^.  = -2°  and  0 = 0°.  This  variation  in  esti- 
mate error  can  be  extended  by  averaging  over  $ and  0^  to  yield  the  average 
interference  - induced  estimate  error  (given  the  presence  of  interference)  vs  off- 
boresight  angle.  Figure  3-3  shows  this  result  for  a cubic  approximation  to 
E(6)  (typical  of  DABSEF)  and  a linear  approximation  for  E(0)  assuming  the  non- 
linear half-angle  processor  of  Eqn.  3. A when  the  signal  to  interference  ratio 
(SIR)  is  10  dB,  and  SNR  = °°. 


The  effects  of  interference  can  be  reduced  somewhat  by  attempting  to 
recognize  pulses  that  were  possibly  interfered  with  and  discarding  them  from 
the  angle  estimation  process.  The  method  employed  for  this  purpose  by  the 
DABS  system  is  to  base  the  reply  azimuth  estimate  only  on  those  pulses 
which  have  a consistent  monopulse  value;  more  detail  is  provided  in 
Appendix  A. 

3.2  ATCRBS  Fruit  Interference 

In  this  case  the  signal-to- interference  ratio  (SIR),  p \ can  vary  over 
the  entire  dynamic  range  of  the  receiver.  From  Figures  like  3-1  and  3-2  one 
can  show  that  the  interference  will  have  a negligible  effect  on  the  bias  error 
if  the  SIR  is  greater  than  +20  dB.  As  the  power  of  interference  increases, 
the  bias  increases,  with  a local  maximum  occurring  at  the  out-of-phase  condi- 
tion. For  SIR's  greater  than  0 dB,  integrating  over  the  phase  difference  leads 
to  a zero  average  value  for  the  monopulse  error.  Since  the  phases  of  the  DABS 
and  ATCRBS  transponders  are  random  from  pulse-to-pulse,  it  is  therefore 
reasonable  to  assume  that  the  greatest  percentage  of  the  sidelobe  fruit  will 
result  in  monopulse  bias  errors  that  will  tend  to  average  to  small  errors. 

For  strong  mainbeam  fruit  or  for  those  few  sidelobe  fruit  whose  power  exeeds 
that  of  the  DABS  transponder,  such  that  the  SIR  is  less  than  0 dB,  the  mono- 
pulse bias  increases  considerably  when  averaged  over  all  phases,  which  indicates 
the  fact  that  the  stronger  interferer  is  starting  to  capture  the  azimuth 
estimator.  When  the  interferer  completely  overpowers  the  DABS  signal,  then  of 
course,  the  monopulse  processor  simply  estimates  the  fruit  azimuth.  Simulation 
and  analytical  studies  have  been  conducted  [7]  which  indicate  that  this  situa- 
tion will  occur  in  only  a small  percentage  of  the  cases  to  be  encountered  in 
the  NAFEC  area  in  1980. 

The  most  important  aspects  of  ATCRBS  interference,  therefore  are  the 
facts  that  when  a fruit  overlap  occurs,  only  a relatively  small  number  of  the 
DABS  bits  will  be  overlapped,  and  of  those  that  are,  since  the  DABS  and  ATCRBS 
transponders  are  incoherent  from  pulse  to  pulse,  the  bias  error  for  low  level 
interference,  (SIR  > 0 dB)  will  scintillate,  and  tend  to  average  out  to  a 
small  value. 

3.3  Specular  Multipath 

In  DABS,  specular  reflections  from  large  flat  surfaces  typically  overlap 
the  direct  received  signal  by  a few  tens  of  wavelengths,  giving  rise  to  errors 
of  the  form  described  by  (3-1).  The  errors  are  constant  for  each  pulse  in  a 
reply  since  the  reflected  signal  is  coherent  with  the  direct  signal. 

For  reflectors  normal  to  the  beam,  the  reflection  azimuth  is  the  same  as 
that  of  the  direct  path  and  leads  to  changes  in  signal  to  noise  ratio,  but 
introduces  no  azimuth  error,  since  Gj.  = 0^,  (in  3-1). 

When  the  surface  is  tilted  slightly  the  effects  of  differences  between 
and  Gg  are  mitigated  somewhat  by  the  reduced  amplitude  of  the  reflected  signal 
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resulting  from  reflection  coefficients  less  than  unity  and  the  cutoff  in  the 
antenna  elevation  pattern-  Only  very  large  and  nearby  structures  will  be 
sources  of  significant  specular  reflection  multipath  entering  the  mainbeam. 
Offboresight  reflectors  producing  sidelobe  multipath  signals  will  have  little 
effect  since  the  relatively  low  sidelobe  gain  helps  to  reduce  errors. 

3.4  Diffuse  Multipath 


The  antithesis  of  specular  multipath  is  completely  diffuse  reflections  in 
which  each  elemental  area  of  the  reflecting  surface  scatters  the  incident  wave 
independently.  Viewed  from  the  interrogator  the  diffuse  reflecting  surface 
consists  of  a large  number  of  point  sources  having  random  amplitudes  and 
phases.  The  diffuse  reflections  combine  randomly  at  the  interrogator  to 
produce  a resultant  that  is  coherent  for  the  duration  of  a transponder  reply 
because  the  change  in  path  length  due  to  a/c  motion  is  insignificant  during 
the  reply.  However,  over  an  ensemble  of  transponder  replies  separated  suf- 
ficiently in  time  or  space,  the  resultant  of  the  diffuse  reflections  has  a 
Rayleigh  amplitude  distribution  and  a random  phase.  Analytical  and  simulation 
studies  of  the  effects  of  diffuse  multipath  on  monopulse  performance  have  led 
to  the  following  observations:  (i)  Time  dispersion  effects  are  not  important 

for  MHz  data  rate  transmission;  (ii)  The  errors  in  azimuth  estimation  are 
negligible  and  orders  of  magnitude  smaller  than  those  that  might  be  observed 
for  specular  multipath  [8].  Therefore,  diffuse  multipath  should  present  no 
serious  limitations  in  direction  finding  capabilities  in  DABS  direction  finding, 
a conclusion  which  has  also  been  supported  by  analysis  of  the  DABSEF  experi- 
mental data. 


3.5  Diffraction  Effects 


Beacon  or  radar  surveillance  systems  measure  target  azimuth  by  estimating 
the  orientation  of  the  wavefront  arriving  at  the  antenna.  DABS  utilizes 
monopulse  processing  to  directly  estimate  the  orientation  of  the  wavefront 
relative  to  antenna  pointing  direction,  while  present  day  ATCRBS  essentially 
computes  the  centroid  of  the  pointing  directions  from  which  the  directional  beam 
receives  signal  energy  greater  than  some  threshold. 

Phenomena  which  perturb  the  wavefront  orientation  must  necessarily  result 
in  errors  in  target  azimuth  estimate.  One  such  phenomena  occurs  when  portions 
of  the  wavefront  are  intercepted  or  blocked  by  obstructions  such  as  buildings. 
The  resulting  field  at  the  ground  receiving  antenna  can  be  shown  to  be  the 
unperturbed  field  less  the  field  through  a rectangular  opening  (with  dimensions 
of  the  blocking  obstacle)  in  an  opaque  screen. 

Recent  investigations  at  Lincoln  Laboratory  (Refs.  10  and  11)  have  shown 
that  in  many  cases  the  dimensions  of  the  obstacle  are  such  that  the  received 
field  can  be  satisfactorily  appoximated  as  the  sum  of  several  rays;  one  is  from 
the  target,  and  the  other (s)  from  appropriate  points  on  the^obstacle  and  having 
appropriate  amplitude  and  phases  relative  to  the  direct  ray  . 

The  effect  of  diffraction  on  present  day  ATCRBS  beamsplitting  can  be 
evaluated  by  modelling  the  azimuth  estimator  as  the  average  of  the  antenna 
pointing  directions  at  which  the  first  and  last  replies  having  a power  greater 
than  some  threshold  are  received.  This  involves  computing  the  received  power 
as  the  antenna  sweeps  past  the  target.  The  azimuth  estimate  error  arises  when 
the  received  power  of  replies  near  the  beam  edge  in  the  absence  of  the  obstacle 
is  modified  by  diffraction  rays  such  that  the  threshold  crossing  azimuths 
change,  leading  to  an  erroneous  azimuth  estimate. 

The  effect  on  DABS  can  be  evaluated  by  superposing  the  various  rays  to 
obtain  the  sum  and  difference  signals  and  then  computing  the  monopulse  pro- 
cessor output.  This  need  only  be  done  when  the  antenna  is  pointing  at  the 
target  since  DABS  makes  angle  estimates  only  near  boresight. 

In  both  DABS  and  ATCRBS  the  azimuth  estimate  error  is  found  to  vary  as  a 
function  of  obstacle  size  and  the  azimuthal  separation  between  the  target  and 
obstacle.  When  the  separation  is  zero,  there  is  no  error,  as  would  be  sugges- 
ted by  symmetry  arguments.  As  the  separation  increases  the  error  oscillates. 

The  oscillation  reflects  the  phase  changes  in  the  obstacle  ray  implied  b the 
differences  between  the  direct  path  to  the  target  and  the  path  through  the 
obstacle.  The  errors  are  zero  when  the  direct  and  diffraction  rays  are  nearly 
in  quadrature. 


* 

Consequently,  one  can  use  simple  interferer  analysis  to  analyze  diffraction 
phenomena . 
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4.0  EFFECTS  OF  ANTENNA  DESIGN 


4. 1 Introduction 


The  antenna  design  parameters  that  have  an  impact  on  monopulse  processing 
performance  are  the  sidelobe  levels  of  the  sum  and  difference  patterns,  the 
vertical  cutoff  of  the  lower  edge  of  the  elevation  pattern, and  the  slope  of  the 
A/E  ratio  [E (0) ] . Since  the  sum  pattern  sidelobes  also  have  an  impact  on  the 
DABS  message  decoding  performance,  it  is  assumed  that  this  is  not  a variable 
parameter.  Therefore,  we  shall  assume  a sum  pattern  which  is  typified  by  the 
standard  28-foot  hogtrough  antenna,  as  shown  in  Fig.  4-1,  and  confine  our 
attention  to  an  examination  of  the  effects  of  difference  pattern  sidelobe 
levels  and  the  elevation  pattern  cutoff.  As  discussed  in  the  preceding  sec- 
tion, the  most  significant  effect  that  interference  has  on  the  monopulse 
performance  is  the  introduction  of  the  bias  term  approximately  given  by  (3-1). 
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The  maximum  error  occurs  at  the  out-of-phase  condition,  <j>  = it  , hence 
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(4-4) 


For  a DABS  target  near  boresight,  Gj  (0g)  = G^O)  and  if  the  sum  and 
difference  patterns  are  normalized  to  G (0) , the  error  expression  can  be  written 
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where  p'  = is  the  reflection  coefficient  after  the  elevation 

pattern  effects  have  been  taken  into  account. 

4 . 2 Mainbeam  Multipath 

For  the  mainbeam  multipath,  0 is  near  boresight  so  G (0  )/G  (0  )wk0 
and  G (0I)  As  1 , so  fl 


(0  - 0 ) < — 

S 1 - p * 


(4-6) 


As  mentioned  in  a previous  section,  the  interference  azimuth  is  related  to  the 
surface  tilt  and  elevation  angle.  The  relationship  is  approximately 


0I  = 0S  + 28ag  (4-7) 

where  8 is  the  surface  tilt  and  the  beacon  elevation,  both  in  radians. 
Therefore,  the  peak  error  is 


~ ~ — pG  (a  )/G  (a  ) 

(0  - 0 ) | < 2- (4-8) 

1 - pGs (ax) /G^ (ag  ) 


Therefore  only  the  elevation  pattern  affects  the  monopulse  performance  in  the 
presence  of  mainbeam  multipath.  However,  its  significance  is  negligible  since 
even  if  we  neglect  the  attenuation  by  this  pattern  and  conservatively  assume  a 
reflection  coefficient  of  0.5,  then  the  peak  bias  error  will  be  .003  for  a 
surface  tilted  at  2 and  a transponder  equipped  aircraft  at  3°  elevation. 

Since  this  error  is  inconsequential,  we  conclude  that  antenna  pattern  design 
is  of  little  importance  in  reducing  the  error  due  to  mainbeam  multipath. 
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Experiment  Results 


4.2.1 

Experiments  were  conducted  at  DABSEF  to  verify  the  above  analytical 
model.  In  these  experiments  the  monopulse  azimuth  error  was  measured  for 
a number  of  radial  test  flights  around  the  DABSEF  facility.  The  standard 
deviation  of  the  monopulse  error  for  a number  of  test  flights  is  summarized 
in  Figure  4-2.  The  data  shown  are  for  two  conditions  of  antenna  elevation 
pointing,  0 elevation  and  -5  elevation.  The  increase  in  monopulse  estima- 
tion error,  when  the  antenna  was  tilted  down  (-5  ) is  likely  due  to  increased 
"in-beam"  multipath. 

4.3  Sidelobe  Multipath 

When  the  specular  multipath  enters  through  the  antenna  pattern  sidelobes, 
then  p'Gj.(0j)  <<  1 and 


(6  - eg)  | = p'  I Gj,  (0I) 
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It  is  clear  that  for  this  case,  the  error  is  smallest  when  the  sum  and  difference 
beam  sidelobes  are  smallest.  When  in  addition  the  target  is  near  boresight, 
the  expression  reduces  to 


lO  -es) 
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(4-10) 


which  shows  that  the  error  can  be  -reduced  by  lowering  the  difference  beam 
sidelobes,  increasing  the  slope  of  normalized  difference  pattern  and  increasing 
the  cutoff  of  the  vertical  elevation  pattern. 

The  direction  finding  (DF)  error  for  a specific  multipath  signal  has  been 
compared  for  two  antennas  of  interest,  namely  the  DABSEF  antenna,  which  has 
independent  E and  A distribution  networks  giving  low  A pattern  sidelobes  as 
shown  in  Figure  4-3,  a,  b,  and  the  RSi  split  hogtrough  shown  in  Figure  4-1. 

The  DF  error  as  a function  of  multipath  signal  azimuth  is  plotted  in  Figure  4- 
4 for  a DABS  target  on  boresight  and  a reflection  coefficient  of  0,5  where  the 
elevation  pattern  attenuation  effects  have  been  ignored.  Figure  4-4  shows  the 
reduction  in  peak  error  resulting  from  the  narrower  mainbeam  of  the 


hogtrough.  (The  DABSEF  3 dB  beamwidth  is  3.2°  while  for  the  RSi  hogtrough  it 
is  2.35  . ) For  larger  multipath  azimuths,  the  peak  error  for  the  hogtrough  is 
larger  because  of  the  higher  difference  pattern  sidelobes  (-14  dB  for  the 
hogtrough,  -30  for  the  DABSEF  antenna)  . 

However,  sidelobe  multipath  errors  are  less  than  peak  mainbeam  multipath 
errors  even  for  the  hogtrough.  When  this  is  considered  together  with  the 
expected  distribution  of  multipath  azimuths  and  amplitudes,  the  usefulness  of 
low  difference  pattern  sidelobes  for  multipath  error  reduction  is  marginal. 

4.4  ATCRBS  Fruit 


The  requirement  for  low  sum  and  difference  pattern  sidelobes  will  be  a 
function  of  sensor  coverage  (terminal  or  en-route) , terrain  shielding,  and 
ATCRBS  fruit  levels.  Curves  showing  the  relation  between  target  aircraft  range 
and  potentially  interfering  ATCRBS  fruit  replies  per  second,  as  a function  of 
antenna  sidelobe  levels,  are  shown  in  Fig.  4-5.  The  fruit  model  used  has 
aircraft  uniformly  distributed  in  range,  total  rate  of  36K  fruit  replies/second, 
and  aircraft  effective  radiated  power  of  +23  dBW. 

Antennas  for  terminal  area  coverage  (50  nmi)  can  have  relatively  high 
sidelobes  (-25  dB  rms) , and  still  not  have  ATCRBS  fruit  replies  affect  mono- 
pulse direction  finding/ (DF)  accuracy.  Simulation  results  indicate  that  16K 
potentially  interfering  fruit  replies/second  will  add  0.07  degrees  of  DF  error 
for  aircraft  interrogated  of f-boresight  at  the  sum/difference  pattern  crossover 
point.  The  terminal  site  is  generally  aided  by  terrain  shielding,  which 
reduces  coverage  and  potentially  interfering  ATCRBS  fruit.  As  can  be  seen  from 
the  curves  shown  in  Fig.  4-5,  en-route  sites  with  good  coverage  will  require 
antennas  with  relatively  low  sidelobes  to  support  accurate  monopulse  DF  for 
aircraft  targets  at  long  range  (to  200  nmi) . 
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5.0  MONOPULSE  PROCESSOR  STABILITY  AND  SENSITIVITY 


The  performance  of  the  monopulse  processor  described  in  Appendix  A has 
been  evaluated  as  a function  of  time,  frequency,  and  amplitude  using  a fixed, 
surveyed,  calibration  transponder.  (See  Appendix  A,  Section  3.3) 

5.1  Stability 

The  stability  data  was  obtained  from  66  calibrations  taken  during  the 
first  6 months  of  1975.  Stability  was  evaluated  at  boresight  and  at  + 1 off 
boresight.  The  results  are  shown  in  Figure  5-1.  The, errors  show  no  long  Q 
term  drift  and  the  standard  deviations  are  .045°,  .021°,  .034  , for  -1  , t).  , 
and  +1  offboresight  respectively. 

5.2  Amplitude  Sensitivity 

Amplitude  sensitivity  was  determined  by  performing  three  calibrations  at 
1090  MHz  while  attenuating  the  received  signal  from  the  calibration  transponder 
by  0 dB,  -23  dB,  and  -39  dB.  The  results  appear  In  Figure  5-2  and  show  that 
the  bias  errors  due  to  amplitude  are  constant  across  the  beam  and  less  than 
0.03°. 


5.3  Frequency  Sensitivity 

The  sensitivity  of  the  monopulse  processor  to  the  received  signal  frequency 
was  evaluated  by  performing  7 calibrations  at  various  frequencies  with  the 
results  shown  in  Figure  5-3.  Note  that  at  1086,  1087,  and  1094  MHz  the  errors 
exceed  0.  1°.  However,  1086  and  1094  MHz  are  outside  the  1090±  3 MHz  National 
Standard  for  ATCRBS  transponders.  Investigations  at  Lincoln  Laboratory,  [9] 
indicate  that  ATCRBS  transponder  frequencies  are  nearly  normally  distributed 
around  1090  MHz  with  a standard  deviation  of  1.75  MHz.  In  Section  3.6  the 
measured  average  azimuth  estimate  error  per  degree  offboresight  for  a large 
sample  of  ATCRBS  targets  is  shown.  Since  the  dominant  error  mechanism  appears 
to  be  transponder  frequency,  Figure  5-4  has  been  drawn  to  depict  the  result 
of  weighting  the  frequency  sensitivity  data  by  the  transponder  frequency 
distribution.  The  expected  average  error  in  Figure  5-4  is  0.02°  per  degree 
offboresight.  The  ATCRBS  azimuth  measurement  is  generally  taken  to  be  the 
average  of  the  estimates  for  replies  that  straddle  boresight,  thus,  reducing 
the  differential  frequency  contribution  to  the  azimuth  estimate  errors  to 
much  less  than  0.1°. 

For  DABS  targets  on  the  roll  call,  the  offboresight  azimuth  at  which  the 
reply  is  received  may  lie  anywhere  within  ±1.5°  (from  scan  to  scan),  but  with 
a tendency  to  lie  nearer  to  -1.5°  (the  leading  edge  of  the  beam)  due  to  the 
nature  of  the  interrogation  scheduling  algorithm.  Under  the  assumption  of  a 
right  triangular  distribution  of  interrogation  offboresight  angles  (favoring 
the  leading  edge),  then  even  a 1087  MHz  target  (having  0.2°/degree  error  slope) 
will  only  experience  a scan-to-scan  standard  deviation  of  0.07°  due  to  this 
effect. 
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6.0  AZIMUTH  ACCURACY 
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The  method  for  evaluating  reply  accuracy  as  a function  of  offboresight  angle 
is  based  on  two  assumptions: 

1.  Aircraft  motion  during  the  beam  dwell  time  is  much 
smaller  than  the  reply  azimuth  estimate  errors. 

2.  An  azimuth  estimate  made  when  the  target  is  at 
boresiglit  would  have  minimum  error. 

The  validity  of  assumption  1 is  guaranteed  be  restricting  the 
analysis  to  targets  at  ranges  greater  than  a few  miles.  Assumption  2. 
is  justified  by  the  frequency  and  amplitude  sensitivities  shown  in 
Section  2. 

Figure  6-1  graphically  illustrates  two  typical  monopulse  characteristics 
where  the  y axis  is  the  monopulse  processor  output  (referred  to  as  the  real 
part),  and  the  x axis  is  the  corresponding  offboresight  angle  in  degrees.  The 
two  tables  differ  in  slope  and  intercept.  Consider  line  A-A  as  the  monopulse 
characteristic  of  the  calibration  transponder  (i.e.,  it  is  the  monopulse  cali- 
bration table)  and  line  B-B  as  the  characteristic  of  an  ATCRBS  target  of  oppor- 
tunity under  surveillance.  This  section  will  explain  the  azimuth  estimation 
errors  induced  by  differences  in  the  two  characteristics  and  describe  the  method 
developed  for  estimating  these  errors  using  reply  data  obtained  during  real 
time  sensor  operation.  (The  errors  shown  in  these  figures  are  exaggerated 
for  clarity.) 


Definitions 
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act 
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est  k 
A0 

act  k 

rP^ 


The  antenna  boresight  azimuth  on  sweep  k. 

The  actual  target  azimuth  on  the  scan. 

The  estimated  azimuth  of  the  reply  on 
sweep  k. 

The  actual  target  of f-boresight  angle  on 
sweep  k. 

The  monopulse  value  (the  "real  part")  of  the  reply  on 
sweep  k. 


A0 


corr  k 


The  calibration  table  value  corresponding 

to  real  part  rp,  , i.e.,  A0  , = C (rp, 
k corr  k.  c * k 

where  C is  the  calibration  table. 


TRAILING  HALF  OF  BEAM 


LEADING  HALF  OF  BEAM 


CORR- 


OFF-BORESIGHT  ANGLE  (DEG) 


CHARACTERISTIC  CALIBRATION  TRANSPONDER 
CHARACTERISTIC  OF  TARGET  TRANSPONCER 

REAL  OFF-BORESIGHT  POSITION  TARGET  WHEN  INTERROGATED 


INDICATED  OFF-BORESIGHT  CORRECTION 


CORR 


6-1.  Indicated  of f-boresight  corrections 


r 


'’l 


Fig.  6.2  shows  the  resultant  azimuth  estimate,  and  individual  reply 
errors^  for  the  case  where  the  target  is  at  10-0  , and  replies  are  received 
at  1.0  intervals  as  the  antenna  scans  from  8.0°  to  12.0  . The  intercept 
of  the  curve  fit  with  the  A0  = 0 Axis,  is  the  reported  aircraft  azimuth,  for  this 
scan  interval.  A normalized  plot  of  reply  errors  vs.  the  estimated  reply 
position  referred  to  boresite,  is  shown  in  Fig.  6-3.  For  the  case  shown,  an 
exaggerated  bias  error  term  was  used,  and  a plot  of  the  true  reply  errors  vs. 
true  position  of  the  replies  referred  to  boresite  is  shown,  for  completeness. 

In  practice  this  bias  term  is  not  measured,  and  is  assumed  small  compared  to 
individual  errors.  Real  data  from  three  scans  are  shown  in  Fig.  6-4,  to 
indicate  the  resultant  output  of  azimuth  estimate  error  as  a function  of 
distance  in  degrees  from  boresite.  A typical  scatter  plot  of  azimuth  estimate 
error  vs.  angle  referred  to  boresite,  is  shown  in  Fig.  6-5,  for  approximately 
100  scans  (500-600  replies) . The  slope  of  the  error  function  for  this  aircraft 
near  boresite  is  approximately  +0.04  /degree.  A measured  distribution  of 
error  slopes  for  approximately  140  aircraft  targets-of -opportunity  is  shown  in 
Fig.  6-6.  For  these  targets,  the  median  error  slope  was  -0.018  /degree,  and 
the  distribution  was  approximately  normal,  with  one  standard  deviation  of 
0 .035°/degree . 
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Fig,  6-2.  Computation  of  estimated  azimuth  errors 


7.0  POSITION  ACCURACY 


7 . 1 Method  of  Deducing  Aircraft  Actual  Trajectory 
Refer  to  Appendix  B. 

7.2  Experimental  Results 
7.2.1  Results  for  ATCRBS  A/C 


Statistics  were  gathered  on  297  altitude  reporting  aircraft  seen  over  6 
days  of  sensor  operation.  A total  of  41,609  reports  were  collected.  The 
overall  standard  deviation  was  .051°,  and  26.4  feet  as  shown  in  Figure  7-1 
and  7-2  respectively.  Figures  7-3  to  7-5  show  the  variation  in  Og  with  ele- 
vation and  azimuth  as  well  as  the  elevation  report  density  function. 

7.2.2  Results  for  DABS  A/C 


Four  DABS  transponders  have  flown  under  DABSEF  Surveillance  for  many 
hours  as  part  of  the  DABS/IPC  (Intermittent  Positive  Control)  testing  program. 
Surveillance  was  usually  maintained  in  both  the  all-call  mode  and  the  discrete 
mode.  The  several  all-call  replies  received  per  A/C  per  scan  were  subjected 
to  the  accuracy  vs  offboresight  angle  analyses  in  the  same  way  ATCRBS  replies 
were.  The  scan-to-scan  sliding  window  was  applied  to  the  position  reports 
based  on  the  single  discrete  reply  per  scan.  The  results  for  about  12000 
scans  of  data  taken  over  several  days  are  shown  below. 
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Figure  7-6  shows  a typical  Intermittent  Positive  Control  (IPC)  encounter 
between  two  DABS  equipped  aircraft.  The  target  reports  are  shown  as  asterisks. 
The  beginning  and  end  of  the  line  segment  associated  with  each  report  represent 
the  smooth  position  and  one  ahead  (4  seconds)  predicted  position  computed  by  the 
IPC  x-y  tracker  (which  assumes  a lo  azimuth  accuracy  of  0.1  degrees).  Thus, 
the  orientation  and  length  represent  tracker  heading  and  speed.  The  high 
quality  of  the  data  can  be  seen  by  referring  to  the  box  of  dimension  d,  cor- 
responding to  three  times  the  0.1  degree  RMS  requirement  for  DABS/IPC  azimuth 
accuracy. 
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Fig.  7-2.  Distribution  of  range  errors. 
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8.0  DIFFRACTION 


A number  of  controlled  experiments  were  conducted  (Refs.  10  and  11)  to 
demonstrate  the  validity  of  the  theoretical  model  for  azimuth  estimation  errors 
caused  by  diffraction  around  shadowing  obstacles.  For  a particular  shadowing 
obstacle  (Hanscom  Field  smokestack,  approximately  10  feet  in  width,  and 
1500  feet  from  the  DABSEF  array  antenna)  we  have  computed  the  expected  azimuth 
estimation  error  for  an  aircraft  at  long  range,  and  low  elevation  angle  (below 
the  top  of  the  tower) . The  error  as  a function  of  the  aircraft  (target)  to 
obstacle  angular  separation  is  shown  in  Fig.  8-1.  The  measured  azimuth  esti- 
mation error  of  a controlled  aircraft  flying  behind  the  obstacle  is  shown  in 
Fig.  8-1  as  well.  It  can  be  seen  that  there  is  extremely  good  agreement  between 
the  expected  and  measured  result.  The  small  differences  are  attributed  to  the 
fact  that  there  are  random  measurement  errors  (approximately  0.04,  one  sigma) 
which  are  convolved  with  the  structured  diffraction  error.  The  theoretical 
model  for  obstacle  shadowing  error  has  been  confirmed  for  difference  geometries, 
using  both  controlled  aircraft,  and  aircraft  targets-of-opportunity . 

The  effect  of  obstacle  shadowing  on  direction  finding  accuracy  can  be 
thought  of  as  the  result  of  distortion  of  the  incident  phase  front  introducing 
error  in  the  estimation  of  direction  of  arrival  of  incident  waves.  Fig.  8-2 
shows  the  theoretical  azimuth  estimation  error  vs.  target-obstacle  angular 
separation,  for  sliding  window  and  monopulse  estimators.  The  magnitude  of 
diffraction  errors  for  typical  shadowing  obstructions  near  airport  beacon 
systems,  is  comparable  to  the  ATCRBS  quantization  error,  and  may  explain  why 
these  structured  errors  have  not  been  observed  in  the  past.  The  structured 
nature  of  errors  induced  by  shadowing  obstructions  is  apparent  in  the  DABS 
data  because  of  the  resolution  of  the  monopulse  estimator.  This  structure 
can  be  seen  from  the  data  shown  in  Fig.  7-5,  where  three  obstructions  are 
clearly  evident.  One  at  agproximately  90  , is  an  FPS-18  radar  antenna 
(700  feet),  another  at  110  , is  a large  semicircular  UHF  array  antenna  (500  feet), 
and  the  third  at  290  , is  the  Hanscom  Field  smokestack  (1500  feet)  referred 
to  above . 
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AZIMUTH  ESTIMATE  ERROR  (DEG) 


9.0  CONCLUSION 


This  report  has  described  the  need  for  a monopulse  system  to  measure 
target  azimuth  using  a single  interrogation/reply  per  scan.  The  requirement 
for  a new  ATC  beacon  antenna  with  an  additional  difference  beam  port  to  pro- 
vide the  monopulse  capability  was  discussed.  Several  methods  of  processing 
the  sum  (I)  omni  (Si)  and  difference  (A)  antenna  outputs  to  produce  an  output 
related  to  target  offboresight  were  reviewed. 

The  choice  of  the  half  angle  processor  because  of  its  unambigous  output 
over  a wide  region  of  the  mainbeam  was  explained.  Various  phenomena  which 
limit  accuracy  or  degrade  performance  were  discussed.  These  included  noise, 
interference  from  reflections,  interference  signals,  and  diffraction  around 
obstructions.  The  effect  of  imperfect  components  in  the  processor  were  con- 
sidered. The  stability  of  the  processor  was  described  along  with  its  sensi- 
tivity to  signal  frequency.  The  sensor  calibration  process  and  its  effects  on 
system  performance  was  described. 

Finally,  data  was  presented  for  a large  number  of  beacon  targets-of- 
opportunity  during  actual  sensor  operation  at  DABSEF.  The  data  showed  reply 
azimuth  estimate  accuracy  as  a function  of  offboresight  angle  and  target  re- 
port accuracy  on  a scan  to  scan  basis. 

A monopulse  direction  finding  capability  has  been  developed,  which  supports 
the  DABS  requirement  (0.1  , one  o)  for  single  reply  direction-of-arrival 
estimation,  a capability  which  also  supports  ATCRBS  direction  finding  at  much 
reduced  PRF  (4-5  replies  per  scan). 
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APPENDIX  A: 


DESCRIPTION  OF  DABS  EXPERIMENTAL 
MONOPULSE  SYSTEM 


1.0  INTRODUCTION 


2.0  DABS  SYSTEM  DESCRIPTION 


2.1  ATCRBS  Mode  of  DABS 


2.1.1 

2.1.2 

2.1.3 

2.1.4 

2.1.5 

2.1.6 


DABSEF  Antenna 

Multi-Channel  Receiver 

Video  Pulse  Quantizer/Digitizer  (VPQ) 

Reply  Processing 

Reply  Correlation 

Report  Azimuth  Determination 


2.2  DABS  Mode 


2.2.1  Reply  Processing 

2.2.2  Position  Measurement  Selection 


A-3.0  MONOPULSE  PROCESSOR 


3. 1 Half- Angle  Implementation 

3.2  Error  Analysis 

3. 3 Monopulse  Calibration 


Figure  A-l. 
Figure  A-2. 
Figure  A-3. 
Figure  A-4. 
Figure  A-5. 
Figure  A-6. 
Figure  A- 7. 


ATCRBS  mode  block  diagram. 

DABSEF  antenna  patterns  vs  azimuth. 

DABSEF  antenna  pattern  vs  elevation. 

DABSEF  antenna  patterns  near  boresight. 
DABSEF  normalized  antenna  pattern. 

Monopulse  processor  error  analysis. 
Monopulse  processor  error  analysis  (Cont.). 
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1.0 


INTRODUCTION 


This  Appendix  contains  an  overview  of  the  DABS  and  ATCRBS  modes  of  DABS  (Section 

2.0)  and  a description  of  the  half-angle  monopulse  processor  (Section  3.0).  The 
material  in  Section  2.0  is  intended  to  provide  the  overall  context  in  which  the 
monopulse  processor  operates,  and  the  various  subsystems  are  described  only  to  a 
level  of  detail  consistent  with  this  purpose. 

2 . 0 _ DABS  SYSTEM  DESCRIPTION 

Section  2-1  will  describe  the  ATCRBS  mode  of  DABS  from  the  antenna  through 
reply  correlation.  (The  final  step,  target-to-track  correlation  is  omitted.) 

Section  2.2  will  describe  the  DABS  mode  starting  with  the  DABS  reply  processor, 
since  the  antenna,  multi-channel  receiver,  and  Video  Pulse  Quantizer  can  be 
considered  common  to  both  modes. 

2.1  ATCRBS  Mode  of  DABS 

Figure  A-l  shows  a block  diagram  of  the  ATCRBS  mode.  Each  of  the  sub- 
systems are  described  in  the  following  sections. 

2.1.1  DABSEF  Antenna 

The  antenna  is  a 512  dipole-element  planar  array  arranged  as  a rectangular  grid 
with  32  columns  of  16  elements  each.  The  elements  are  printed  circuit  dipoles. 

Each  vertical  column  of  dipoles  is  housed  within  its  own  radome.  The  antenna 
support  structure  is  a lightweight  aluminum  frame  designed  for  low  azimuth 
sidelobes.  The  antenna  also  provides  separability  of  networks  for  elevation  and 
azimuth  distribution  (and  individual  sum  and  difference  networks).  Table  A-l 
summarizes  the  significant  antenna  performance  characteristics. 

These  include: 

Operating  Frequency  Bandwidth  - Largely  because  of  the  wideband  design  approach 
adopted,  the  results  show  the  antenna  to  have  considerable  operating  bandwidth 
with  sharp  elevation  pattern  cut-off  and  low  sidelobes  over  the  full  band. 

Low  Azimuth  Pattern  Sidelobes  - Over  the  range  of  +1°  to  +35°  in  elevation,  the 
maximum  level  of  all  azimuth  sidelobes  measured  for  both  sum  and  difference 
patterns  (0°  to  360°  in  azimuth)  was  26  dB. 

Elevation  Coverage  Sector  - Coverage  over  35°  with  uniformity  of  level  to  within 
±1  dB  was  achieved. 

Monopulse  Characteristics  - Difference  pattern  with  highly  symmetrical  charac- 
teristics and  null  depth  of  greater  than  40  dB. 


TABLE  A-l  DABSEF  ANTENNA  CHARACTERISTICS 
Array  Antenna  for  DABS  Experimental  Facility 
Main  Array  Antenna 


Frequency  of  Operation 

Aperture  Size 

Number  of  Array  Elements 

Sum  and  Difference  Patterns 
for  Monopulse 

Azimuth  Sum  Beamwidth 

Azimuth  Difference  Null  Depth 

Azimuth  Pattern  Sidelobes  (5°  EL) 

Elevation  pattern  with  sharp 
cutoff  at  horizon 

Elevation  pattern  slope  at 
-6  dB  point 

Elevation  Sector  Coverage 

Antenna  Input  VSWR 
RF  Power-handling  Capability 

Gain  - sum 

- difference 
Polarization 
Cross-Polar ization 
RF  Losses 


970  to  1150  MHz 
8'  by  20' 

512  dipoles 


3 0 

40  dB  max 

-30  dB  max  at  1090  MHz 

3. 1 dB/degree  at  1090  MHz 

Uniform  within  2 dB  over 
35° 

1.4:1 

5 kw  peak 
1 kw  avg. 

22.5  dB 

20.5  dB 
Vertical 

-40  dB 
1.7  dB 


Omni  Antenna 


Elevation  Pattern 
Azimuth  Pattern 


Matched  to  Array  antenna 

Uniform  around  360°  to 
within  3 dB  at  +5  elevation 
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Omni-Antenna  Performance  - Elevation  pattern  same  as  array  pattern  and  azimuth 
coverage  uniform  to  within  ±1.7  dB. 

Figures  A-2  through  A-4  show  the  azimuth  sum  and  difference  patterns,  elevation 
pattern,  and  sum  and  difference  near  boresight.  The  three  outputs,  sum  (E), 
difference  (A),  and  omni  (ft) , are  processed  by  the  multichannel  receiver. 

2.1.2  Multi-Channel  Receiver 

The  receiving  system  consists  of  a front  end  and  IF/video  processor.  The 
front  end,  located  near  the  antenna,  converts  £ and  A to  E+jA  and  A+jE  (A 
and  B respectively)  limits  A,  B,  and  the  omni  (ft),  converts  to  IF,  amplifies, 
and  sends  them  to  the  IF/vodeo  processor.  The  IF/video  processor  filters  A, 

B,  and  ft  and  forms  the  monopulse  video  and  log  |e|  video  from  A and  B,  and 
forms  log  |ft|  video.  Finally  the  video  quantizer  operates  on  log  |e|,  log 
|ft|,  and  log  J A J to  yield: 

a)  QEA  which  indicates  when  log  |e|  video  exceeds  a threshold 

b)  QEPS  which  indicates  when  log  J E | video  has  a positive  slope 

c)  QENS  which  indicates  when  log  (e|  video  has  a negative  slope 

d)  QSLS  which  indicates  when  log  j E | video  exceeds  log  |ft|  or  log 

| A | video  by  a certain  amount.  This  signal  represents  the  Receive 
Sidelobe  Suppression  (RSLS)  function  which  restricts  processing  to 
the  mainbeam. 

Outputs  (a)  to  (d)  plus  the  monopulse  video  (described  in  Section  3.0)  are 
input  to  the  video  pulse  quantizer  digitizer. 


2.  1.3  Video  Pulse  Quantizer /Digitizer  (VPQ) 


The  VPQ  samples  QEA,  QEPS,  and  QENS  at  8.276  MHz  and  logically  operates 
on  the  samples  so  as  to  recognize  the  occurrence  of  pulse  edges.  The  VPQ 
outputs  a data  stream,  QLE,  clocked  at  8.276  MHz,  that  has  a "1"  whenever  a 
pulse  leading  edge  is  detected  and  "0"s  otherwise.  Leading  edges  are  declared 
whenever  log  | E | video  exceeds  a threshold  (QEA  is  "on")  and  the  slope  of  log 
| E | video  transitions  from  high  to  low.  This  logic  enables  declaration  of 
leading  edges  for  overlapped  pulses  when  the  second  pulse  is  sufficiently  high 
in  amplitude.  Trailing  edges  are  recognized  in  a similar  fashion.  They  are 
used  to  reject  short  pulses  and  declare  leading  edges  in  certain  pulse  overlap 
conditions.  The  monopulse  video  is  sampled  and  quantized  to  8 bits  shortly 
after  each  leading  edge  (except  those  in  overlap  conditions).  The  QSLS  signal 
is  also  sampled  after  each  leading  edge  resulting  in  the  SLSF  clocked  data 
stream  (also  8.276  MHz). 
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2.1.4  Reply  Processing 

The  DABSEF  ATCRBS  reply  processor  searches  the  ZLE  data  stream  for  pulses 
separated  by  the  Fj,  F2  spacing  (requiring  that  at  least  one  is  main  beam).  The 
bracket  detection  logic  is  capable  of  resolving  various  types  of  interreply 
garbling  and  interference.  When  F^jF^  brackets  are  detected  the  associated  code 
pulse  positions  in  the  2LE  stream  are  examined.  The  presence  or  absence  of  a 
leading  edge  in  the  vicinity  of  those  positions  results  in  a "1"  or  "0"  being 
declared.  Due  to  noise,  garbling,  and  other  phenomena,  it  is  possible  for  these 
declarations  to  be  in  error.  Therefore,  to  aid  in  subsequent  processing,  each 
code  bit  has  a corresponding  confidence  bit.  The  setting  of  these  bits  depends 
on  the  consistency  of  the  monopulse  values  for  the  leading  edges  and  the  presence 
or  absence  of  edges  in  the  sample  positions  adjacent  to  that  of  the  nominal  code 
pulse  positions.  The  monopulse  value  of  each  pulse  is  compared  to  the  "average" 
of  the  preceding  pulses.  The  average  is  initialized  with  the  F^  value  (except 
during  garbling  when  F?  is  used).  It  is  updated  by  adding  the  current  average 
to  the  new  pulse's  value  and  dividing  by  two,  if  the  new  pulse  is  sufficiently 
close  to  average.  If  not,  the  old  average  is  retained  and  the  code  bit  is 
declared  a low  confidence  "1".  When  decoding  and  monopulse  averaging  is  com- 
plete, the  reply  is  sent  to  the  reply  correlation  function  wi  .1 


a) 

mode  indicator 

b) 

range 

c) 

monopulse  average 

d) 

antenna  boresight  at  bracket  detection  time 

e) 

code 

f) 

code  confidence 

2.1.5 

Reply  Correlation 

Replies  from  successive  sweeps  are  correlated  on  the  basis  of  range, 
azimuth,  and  code  and  combined  into  target  reports.  A reply  that  does  not 
correlate  with  replies  from  previous  sweeps  is  considered  an  initial  target 
report  having  the  range,  azimuth,  and  code  (Mode  A or  C)  of  the  reply.  The 
rules  for  updating  target  reports  with  correlating  replies  result  in  the  report 
range  and  azimuth  being  the  average  of  the  two  replies  that  straddle  boresight 
or,  if  all  replies  are  from  the  same  side  of  boresight,  the  values  for  the  reply 
closest  to  boresight.  The  report  codes  (A  and  (!)  and  code  confidence  words  are 
obtained  by  logical  operations  on  the  reply  codes  and  code  confidence/ 
bits . 
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2 . 1 . 6 Report  Azimuth  Determination 


The  averaging  of  sampled  monopulse  values  during  reply  processing,  combined 
with  the  reply  averaging,  reduce  the  noise  variance  in  the  ATCRBS  target  re- 
port. 

In  reply  processing  the  "add  and  divide  by  two"  process  yields  a sigma  due 
to  noise  of 


reply  ^ pulse51 


3 1 + 


This  is  a weak  function  of  n,  and  since  n can  be  as  small  as  two  (brackets 
only),  ' will  be  taken  as  apulse/\2.  The  averaging  of  the  two  reply  azimuths 
that  straddle  boresight  in  reply  correlation  provides  another  factor  l/v2,  so 


report 


(for  noise  effects) 


2.2  DABS  Mode 


2.2.1  Replv  Processing 


The  DABS  reply  processor  operates  on  the  VPQ  outputs  to  generate  a 16  MHz 
data  stream  (,SQ^D)  indicating  the  presence  of  pulses  above  the  DABS  threshold. 

SQ  D is  then  input  to  the  DABS  preamble  detector  which  searches  for  preamble 
pulse  sequences.  When  a DABS  reply  preamble  is  detected  a trigger  pulse  is 
generated  which  causes  SQZD,  QSLS,  and  log  Jz  to  be  sampled.  The  sampling 
rate  is  twice  the  information  bit  rate  to  enable  decoding  of  the  two  Pulse 
Position  Modulation  "chips"  that  make  up  each  information  bit.  The  monopulse 
video  is  sampled  in  the  same  way  and  converted  to  8 bit  digitized  values. 

Finnallv,  the  log  Z is  sampled  once  per  chip  (short  reply.'  56  bits,  112  chips). 
The  two  log  chip  samples  per  information  bit  are  then  compared,  and  a data 

stream  at  the  information  bit  rate  is  output  indicating  which  of  the  two  chips  had 
a greater  amplitude. 

The  -essage  bit  processor  then  logical lv  operates  on  the  log  jzj  , OLSL  and 
the  -onopulse  video  chip  samples  plus  the  amplitude  comparisons  to  declare  the 
binary  value  of  each  information  bit.  As  with  ATCRBS,  a confidence  value 
jc  omnanies  each  declaration.  The  message  bits  are  then  error  detected  and 

corrected . 

2.2.2  P sit  ion  Measurement  Selection 

DABS  nos  it  ion  measurements  are  made  on  the  basis  of  only  one  reply  per  scan 
ever  if  others  were  received.  The  general  rult  is  to  use  the  roll-call  reply 
closest  to  boresight  or,  In  the  absence  of  roll-call  replies,  to  use  the  all- 
: ; ; renlv  having  the  shortest  range. 
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3.0  Monopulse  Processor 

3.1  Half  Angle  Implementation 


The  monopulse  processor  implemented  at  DABSEF  is  the  half  angle  realization 
described  in  Section  2.3  of  the  Report.  Two  independent  measurements  of  tan-^  (A/Z) 
are  made  and  summed  so  that  the  processor  output  is  2 tan  A/Z  (for  a triangular 
detector) . 

Since, 

2 tan  ' (A/Z)  = [arg  (Z  + jA)  - arg  (Z  - jA)]  = 0 
the  optimum  Bell  Labs  result  is  achieved. 

-i 

The  measurements  of  tan  (A/Z)  are  obtained  by  phase  detecting  between  (Z  + jA) 
and  jZ  in  the  one  channel  and  between  -(A  + jZ)  and  -Z  in  the  other.  This 
eliminates  the  ambiguities  incurred  when  0 exceeds  90  (i.e.,  A = Z)  in  the 

Bell  Labs  approach.  Utilizing  both  channels  avoids  any  loss  in  signal  to  noise 
ratio,  and  reduces  the  effect  of  phase  errors  in  the  hybrid  channels  by  one 
half. 


3.2  Error  Analysis 

This  section  computes  the  azimuth  estimation  accuracy  of  the  DABS  half 
angle  processor  and  the  DABSEF  antenna  described  in  Section  2.3.  We  age  concerned 
with  a single  DABS  reply  within  the  antenna  half  power  beamwidth  (±1.5  ),  and 
in  an  elevation  angle  to  A/C  regime  limited  to  the  coverage  region  0 to  5 . 

In  this  region  we  assume  that  the  monopulse  pattern  is  constant. 

The  A/Z  ratio,  of  the  DABSEF  antenna,  vs  angle  referred  to  boresight,  is 
shown  in  Fig.  A-5a.  The  derivative  of  this  curve,  in  A/Z  units  per  azimuth 
degree,  is  shown  in  Fig.  A-5b.  Tnese  curves  are  used  to  convert  processor 
errors  in  A/Z  units  to  azimuth  errors  in  degrees  vs  angle  referred  to  boresight. 

The  sources  that  contribute  to  errors  in  estimated  angle  of  arrival  consist 
of,  but  are  not  limited  to,  the  following: 

1.  Differential  pre-combiner  gain/phase. 

2.  Differential  channel  phase  error  due  to  transponder  frequency 
range  (±3  MHz),  and  dynamic  range  of  received  signals  (-20 
dBm  to  -79  dBm) . 

3.  Phase  error  due  to  sampling  the  pulsed  transient,  processor 
output  signal,  in  a non-quiescent  state. 

4.  Noise. 
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DABSEF  ANTENNA  A/I  RATIO 
(a) 


0,0  1.0  2,0  3,0 

OFF-BORESIGHT  ANGLE  (DEG) 


DABSEF  ANTENNA  SLOPE  OF  A/I  RATIO 
(b) 


Fig.  A-5.  DABSEF  normalized  antenna  pattern. 


Other  possible  sources  of  estimation  error  relate  to  the  antenna  system, 
coupling  of  shaft  encoders  to  the  driven  system,  mechanical  stability,  quanti- 
zation error  in  reading  antenna  position,  and  quantization  error  in  recording 
the  real  part  of  the  measured  A /E  ratio.  These  are  considered  to  be  small 
effects,  and  are  not  considered  here.  Another  source  of  estimation  error  relates 
to  the  change  in  the  monopulse  channel  levels  vs  A/C  (target)  elevation  angle 
(approximately  Cos  EL  dependent),  but  is  not  considered  here. 

In  each  case  a standard  deviation  for  the  perturbing  influence  is  assumed 
and  transformed  to  the  standard  deviation  in  the  A/S  ratio  by  making  usa  of  the 
relationship  between  the  processor  output  and  the  perturbation.  Let: 


enm’D d = Perturbing  influence 
r K 1 Kd 


ePFTRB  mean  °f  SPRTRB 


ePRTRB  = standard  deviation  of  epRTRB 


Then : 


A/E  = a given  A/E  input  ratio 


IN 

A/Equi  = processor  output  A/E  ratio 

f (...)  = the  function  A/Equt  = f (A/EIN,  epRTRB) 

o , output  standard  deviation  (a  function  of  A/E_„  and 

A/E n - IN 

PRTRB,  and  a ). 

ePRTRB 


’*V  “ %rtrb  *nmi> 


a ^ /„  is  easily  converted  to  off  boresight  angle  estimation  error  for  a 

given  of f^&oresight  angle  by  utilizing  the  antenna  A/E  ratio  versus  off  boresight 
characteristic . 


Let : 


9 = off  boresight  angle 
A9  = off  boresight  angle  estimation  error 
0^(0)  = standard  deviation  of  A0  for  a given  9 

g(.)  = the  function  A/EpN  = g(9) 
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Then: 


1_ 

,n,  ^ d ( g ( 0 ) ) 
°A0(0)  = ~ 3T~ 


A/Iout  (g(9» 


or  Combining  (1)  and  (2) 

_1 

o (9)  = dg(B)  f'(g(r),  e ) a 
A9  dT  PRTRB  ePRTRB 

Starting  at  the  front  end  of  the  processor,  the  first  possible  source  of 
error  is  differential  gain/phase  variations  between  the  A and  E channels,  ahead 
of  the  tt/2  combining  hyl  id. 


Let: 


g = ratio  of  A channel  to  E channel  gains  ( 3 * in  dB) 
g = 1 (0  dB) 
oOI=  .25  dB 


Then 


f (A/E , g)  = B A/E 
a 


a/e 


OUT 


oB  4 n 


Since: 

g 


rat  Lo 


= 10 


(B'/20) 


Then: 

0g  . 

ratio 

Therefore 

a 


10  (ag  ^20)-l 


A/E  = f (10(',p,/20)  -1) 


= (.0202)  A/E 

A plot  of  this  error  function  is  shown  in  Fig.  A-6a,  where  scales  have  been 
transformed  to  antenna  units,  using  Figs.  A-5a,  A-5b. 

Let  <j>  = phase  difference  between  the  A and  E channels 

^ = 0 


Then  : 


f (A/E  , 4>)  = 


1 


tan  tan 


- 1 A cos  4' 

E - A sin  <}i 


+ tan  tan 


-1  A cos 

E + A sin 


The  derivative  evaluated  at  <j>  = 0 is  zero.  Consequently  the  processor  is 
insensitive  to  this  source  of  error. 
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e.  DUE  TO:  c AND  d 

Fig.  A-6.  Monopulse  processor  error  analysis. 
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System  phase  errors  consist  of  errors  which  occur  principally  due  to  mis- 
matched channel  $ (w)  transfer  functions,  and  those  which  are  introduced  by 
active  elements  as  a result  of  change  in  with  signal  level.  The  elements 

of  the  processor  which  are  in  the  first  category  are:  the  front-end  pre-select 

band  pass  filters,  differential  channel  line  length  (Ax),  the  IF  band  pass 
filters,  and  the  A/E  sampling  A/D  converter  (transient  phase  response  varies  with 
input  frequency  due  principally  to  the  IF  band  pass  filters). 

The  output  of  the  phase  detectors  would  be  exactly  one  half  the  angle 
between  2 + jA  and  A + jE  for  ideal  triangular  detectors.  Thus,  the  detector 
output  for  a phase  error  <f>  is 

tan_1(A/SIN)  + | 

The  corresponding  indicated  a/£qjjX  is  then 

A/I0UT  = tan  (tanl  A/5:IN  + 2} 

Differentiating  and  evaluating  at  <j>  = 0 yields 

A/zoux  = £i  (1  + (A/E)2) 

Peak-to-peak  phase  errors  for  the  causes  cited  will  be  less  than  8°.  It 
is  assumed  that  the  mean  error  is  zero  (at  1090  MHz),  and  that  transponder 
frequencies  will  be  normally  distributed  about  1090  MHz.  A la  phase  error  of  2 
is  chosen,  a plot  of  the  error  function  is  shown  in  Fig.  A-6b , where  the  scales 
have  been  transformed,  as  before.  In  practice,  the  effect  of  imperfect  tracking 
of  the  IF  band  pass  filters  with  signal  frequency  can  be  minimized  by  careful 
selection  of  the  processor  output  sampling  time  referenced  to  the  pulse  leading 
edge. 


The  IF  constant  phase-limiters  are  the  only  elements  in  the  processor 
which  introduce  phase  error  as  a function  of  signal  level.  The  limiters  in  the 
(E  + jA)  and  the  (A  + jE)  channels  do  not  track  in  phase  over  the  required  60  dB 
dynamic  range.  This  is  a minor  effect,  and  peak-to-peak  tracking  phase  error  is 
expected  to  be  less  than  3 . It  is  "assumed"  that  received  signal  levels  will 
be  normally  distributed  about  a median  (calibration)  level,  and  the  la  phase 
error  is  therefore  taken  to  be  0.75  . The  error  in  the  A/E  ratio  is  then 

cA/2 (Limiter  Phase)  = .0065  (1  + |a/e|2)  (5) 

A plot  of  this  error  function  is  shown  in  Fig.  A-6c. 

The  limiter  output  level  will  vary  as  a function  of  input  level.  It  is  expected 
that  peak-to-peak  variation  will  be  less  than  0.75  dB  over  the  60  dB  input 
signal  range.  Again,  assuming  normally  distributed  signal  levels  about  a 
median  level,  a la  output  level  variation  of  approximately  (.75/4)  dB  has  been 
-chosen. 
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The  detector  output  for  a differential  gain  ratio  of  8 between  £ + jA  and 
A + jl  channel  is 


tan[8  tan  (A/I)] 

Differentiating,  evaluating  at  6 = 1 and  expressing  o 
yields  ratio 


in  terms  of  ar 


°a/£  = tan_1  (A/E)  -^rio  • [1°(Y/20)  -n  [1  + lA/si2] 

Y = la  differential  output  level  in  dB 
= .75  dB/4  (estimate) 

/ _ 1 2 

0 . (Limiter  Output  Level)  = 3.81  x 10  tan  (A/I)  . (1  + |A/l]  ) 


A plot  of  this  error  function  is  shown  in  Figure  A-6d,  and  the  combined  error 
function  (r.s.s.),  is  shown  in  Figure  A-6e. 

The  last  term  which  affects  the  estimate  of  A/I  is  additive,  white, 
Gaussian  noise.  The  la  error  in  A/I  due  to  noise  (SNR)  per  pulse  is, 


aA/i (Noise) 


(1  + | A/I | ) 


where  — z — r-  accounts  for  variations  in  SNR  across  the  beam  due  to  the 

1(8) 

sum  beam  pattern.  Plots  of  this  function  are  shown  in  Figure  A-7a  for  +30 
and  +42  dB  SNR  with  scales  transformed  as  before.  Curve  No.  1 was  combined 
with  other  error  sources  in  a r.s.s.  sense  to  obtain  the  curve  shown  in 
Figure  A- 7b.  This  represents  the  expected  estimation  error  performance  of 
the  processor  and  DABSEF  antenna  for  DABS  replies  at  or  above  +42  dB  referred 
to  MUSL. 


3.3  Monopulse  Calibration 

The  DABS  sensor  is  calibrated  against  a fixed  transponder (s)  at  a known 
azimuth (s) . The  result  is  a look-up  table  which  relates  digitized  monopulse 
video  samples  (y)  to  offboresight  angle.  This  section  will  describe  the 
calibration  process. 
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Fig.  A-7.  Monopulse  processor  error  analysis  (Cont.). 


Let: 


— — 


0 = Azimuth  of  target  (or  calibration  source  during 

formation  of  look-up  table)  with  respect  to  some 
reference,  for  example  true  north. 

BS  = azimuth  of  difference  beam  null  with  respect  to 

north.  (The  difference  beam  null  is  an  arbitrary 
but  convenient  definition  of  beam  azimuth). 

bs.  , = The  output  of  the  antenna  shaft  encoder 

xnd 

A0  = The  angle  between  the  target  and  the  null  of  the 

difference  beam. 


Y = digitized  half-angle  processor  output 


Then: 

0 = BS  - A0 

true  true  true 


(1) 


Further  let: 

A0 

true 


= f (y) 


(2) 


and  define  f(.)  for  an  "aligned"  system  where  alignment  corresponds  to  a mono- 
pulse bipolar  video  output  of  zero  (corresponding  to  y = 128)  when  the  target 
is  at  the  difference  beam  null.  The  function  f(.)  is  then  the  equivalent  of 
Equation  2-6  of  Section  2.3. 

Now  introduce  two  fixed  offsets,  one  in  the  A/D  correcter  that  digitizes 
the  monopulse  video  and  the  other  in  the  azimuth  shaft  encoder. 


Then 


6 = A/D  offset 


ABS  = bs.  - BS 

ind  true 


Y = f 1 (A0  ) + 6 

true 


so  that  combining  (1),  (3),  and  (4) 


0 = bs.  - (f(y-6)  + ABS) 

true  md 


(3) 

(A) 

(5) 
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where  the  term  in  parenthesis  is  obviously  the  monopulse  look  up  table  (a  func- 
tion of  y)  which  is  interpreted  as  indicated  off  boresight  angle  (A6),na 


0 = bs  . , - A0 . , 

true  ind  ind 


A9ind  = C(>>  = f + ABS 


(6) 

(7) 


Finally,  f(y-<5)  is  the  angle  between  the  target  and  the  difference  beam 
null  and  ABS  is  the  angle  between  the  difference  beam  null  and  the  processor/A/D 
system  null. 

The  look-up  table  C(y)  is  formed  using  replies  from  a surveyed  target; 
i.e.,  0 is  known.  Each  reply  (i)  contains  a value  bs . and  y.,  which  are 

combined  as  shown  below  to  obtain  C(y).  1 


C(Y)  = 


E 


(bs.  - 0 ) g (y  - y) 

i l true  6 ' i ' ’ 


E 


g (y  • - y) 


where  g(Y|  - y)  = 0 Y^  ^ Y 


Y.  = Y 


Note  that  the  offsets  6 and  ABS  displace  the  look  table  along  the  y and  A0 
axes.  However,  as  long  as  the  offsets  are  unchanged  from  their  values  at 
calibration  time,  the  azimuth  estimates  are  unbiased. 

If,  however,  the  survey  of  the  calibration  source  is  in  error  a bias  is 
introduced  into  all  position  measurements.  This  bias  can  affect  ATC  system 
configurations  that  use  data  from  more  than  one  sensor. 

For  any  given  target  at  any  given  time  and  position  there  exists  a cor- 
responding lookup  table  C (y)  which  relates  y to  A0  . CT(y)  takes  into 
account  all  target  characteristics  (such  as  transporter  frequency  and  elevation) 
as  well  as  the  instantaneous  characteristics  of  the  antenna,  receivers,  half 
angle  processor,  etc.  The  consequence  of  using  C(y),  the  stored  lookup  table, 
rather  than  CT(y)  is  to  create  an  error  in  the  azimuth  estimates 
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APPENDIX  B:  METHOD  OF  DEDUCING  AIRCRAFT 

ACTUAL  TRAJECTORY 

A fundamental  objective  of  many  data  analysis  programs  is  estimation 
of  an  actual  aircraft  trajectory,  given  a sequence  of  nosition  measurements 
for  that  aircraft.  The  sliding  window  least  squares  curve  fit  is  one  tech- 
nique for  achieving  such  an  estimate.  The  fundamental  characteristic  of  the 
procedure  is  the  use  of  a sliding  window  to  obtain  a position  estimate  cor- 
responding to  a target  report.  The  window  defines  all  data  points  contained 
in  a fixed  time  period  enclosing  the  target  report  for  which  an  estimate  is 
desired.  In  this  application  a second  order  least  squares  curve  fit  is  applied 
to  all  applicable  data  points,  and  the  "true"  position  is  calculated  from  the 
curve  fit  equations.  As  the  procedure  advances  to  the  next  target  report,  the 
time  window  is  also  moved.  This  produces  a sliding  window  which  defines  what 
data  is  to  be  used  in  the  curve  fit  procedure. 

To  illustrate  the  sliding  window  procedure,  assume  an  estimate  of  the 
target  position  is  desired  for  scan  N.  Target  report  data  (range,  azimuth  and 
time)  for  scans  N-5  to  N+5  are  used  in  the  fit,  (assuming  the  time  window  is 
selected  to  correspond  to  11  scans).  A least  squares  curve  fit  is  determined 
from  the  target  reports,  yielding  equations  for  range  and  azimuth  as  a function 
of  time: 


p(t)  - al  + a2  (t  - tN)  + a3  (t  - tn)2 
0(t)  = bl  + b2  (t  ' tN)  b3  (t  " tn) 2 

In  the  equations,  t denotes  the  time  associated  with  the  target  report  for 
scan  N,  and  the  least  squares  estimate  of  the  aircraft  range  and  azimuth  would 
be  a^  and  b.  for  scan  N.  Velocity  and  acceleration  components  are  also  avail- 
able if  heading  and  ground  speed  estimates  are  desired. 

In  an  actual  implementation,  it  is  not  guaranteed  that  target  reports 
exist  for  each  scan.  If  the  number  of  data  points  in  the  window  is  less  than 
some  minimum  (typically  7),  a fit  is  not  attempted.  In  addition  to  missing 
target  reports,  insufficient  data  may  also  exist  if  target  report  measure- 
ments are  flagged  as  "bad"  data  not  to  be  used  in  the  curve  fit.  However,  if 
the  measurements  for  a given  scan  are  rejected  and  sufficient  data  is  still 
available,  a curve  fit  is  still  executed. 

The  bad  point  rejection  scheme  is  applied  independently  to  each  curve  fit 
(range  and  azimuth),  so  partial  position  estimates  may  exist  for  some  scans. 

If  the  actual  aircraft  trajectory  was  the  desired  program  output  this  would  not 
be  useful,  but  in  this  application  a partial  position  estimate  can  still  be 
used  to  calculate  a measurement  error  in  one  of  the  two  measurement  coordinates. 
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The  use  of  independent  range  and  azimuth  curve  fits  is  designed  to  maxi- 
mize the  amount  of  measurement  error  information  achievable  from  the  input 
data.  Independent  curve  fits  in  range  and  azimuth  are  also  desirable  in  the 
sense  that  poor  data  in  one  coordinate  does  not  corrupt  the  quality  of  the 
curve  fit  for  the  other  coordinate.  For  a "typical"  aircraft,  the  standard 
deviation  of  the  range  measurement  error  is  20  -25  feet,  while  the  standard 
deviation  of  the  azimuth  measurement  error  is  .04°.  At  a range  of  35  miles, 
the  variance  of  the  azimuth  measurements  would  be  about  150  feet,  so  it  is 
clear  that  azimuth  errors  are  the  dominant  factor  in  the  surveillance  errors. 

If  an  X - Y curve  fit  were  used,  azimuth  errors  would  translate  into  the  X and 
Y curve  fits,  and  could  produce  apparent  range  error  while  simultaneously 
reducing  the  azimuth  errors.  The  coupling  that  results  from  an  X-Y  curve  fit 
was  the  primary  factor  in  selecting  the  range-azimuth  fit  for  the  programs. 

In  contrast,  it  is  noted  that  the  second  order  equations  in  range  and 
azimuth  cannot  model  the  idealized  cases  normally  assumed  in  similar  curve  fit 
programs.  If  an  aircraft  were  flying  with  constant  acceleration  on  a perfectly 
linear  track,  ideal  data  applied  to  an  X-Y  curve  fit  would  produce  no  error. 

The  same  data  applied  to  a range-azimuth  curve  fit  produces  incorrect  indi- 
cations of  errors  in  the  measurements.  Comparisons  of  the  two  tehchniques  in 
simulations  and  with  actual  data  indicate  that  the  differences  in  the  two 
techniques  result  in  relatively  small  differences  when  averaged  over  many 
scans.  In  worst  case  geometries,  the  differences  were  typically  less  than  10 
feet  in  range  and  .02  in  azimuth  for  aircraft  at  ranges  greater  than  10  miles 
from  the  sensor  operating  at  moderate  speed  (less  than  200  knots)  on  basically 
linear  trajectories. 

In  selecting  a range-azimuth  curve  fit,  it  was  concluded  that  the  accuracy 
was  sufficient  for  estimating  the  cumulative  measurement  errors.  It  also  lends 
itself  better  to  bad  point  rejection  schemes  since  the  entire  measurement  does 
not  have  to  be  rejected  if  the  azimuth,  for  example,  appears  inconsistent. 

Finally,  a turn  detection  in  X-Y  is  incorporated  and  used  to  filter  out 
curve  fit  residuals  for  scans  on  which  the  aircraft  is  turning. 

Each  residual  is  incorporated  into  two  data  bases. 

1.  It  is  combined  with  other  residuals  from  the  aircraft. 

The  output  is  the  standard  deviation  of  azimuth  estimate 
error  for  that  aircraft. 

2.  It  is  combined  with  the  residuals  for  all  other  aircraft. 

The  output  is  the  standard  deviation  of  all  the  residuals, 
and  standard  deviation  plotted  vs  azimuth  and  vs  elevation. 
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